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Neutrinoless double-/? decay is of fundamental importance for determining the neutrino mass. 

Although double electron [I3~P~) decay is the most promising mode, in very recent years interest 
in double positron (/?''"/?^) decay, positron emitting electron capture {EC/3'^), and double electron 
capture (ECEC) has been renewed. We present here results of a calculation of nuclear matrix 
elements for neutrinoless double-/?"*' decay and positron emitting electron capture within the frame¬ 
work of the microscopic interacting boson model (IBM-2) for ®®Ni, ®'*Zn, ^®Kr, ®®Ru, *^°®Cd, 

and *^®®Ce decay. By combining these with a calculation of phase space factors we calculate 
expected half-lives. 


PACS numbers: 23.40.Hc,21.60.Fw,27.50.+e,27.60.+j 

I. INTRODUCTION 

Double-/? decay is a process in which a nucleus {A, Z) 
decays to a nucleus {A, Z ± 2) by emitting two electrons 
or positrons and, usually, other light particles 

{A, Z) —>■ {A, Z ±2) 26^^ anything. (1) 

Double-/? decay can be classified in various modes ac¬ 
cording to the various types of particles emitted in the 
decay. The processes where two neutrinos are emitted 
are predicted by the standard model, and 2vj3~ decay 
has been observed in several nuclei. For processes not al¬ 
lowed by the standard model, i.e. the neutrinoless modes: 
Oz//?/?, Qv(3EC, OvECEC, the half-life can be factorized 
as 

^ = Go.\Mo,\^\f{m,,Uei)f, (2) 

where Gqii is a phase space factor, Mq^, is the nuclear ma¬ 
trix element, and firrii, Uei) contains physics beyond the 
standard model through the masses and mixing ma¬ 
trix elements Uei of neutrino species. For all processes, 
two crucial ingredients are the phase space factors (PSFs) 
and the nuclear matrix elements (NMEs). Recently, we 
have initiated a program for the evaluation of both quan¬ 
tities and presented results for /?“/?“ decay (iH3. This 
is the most promising mode for the possible detection of 
neutrinoless double-/? decay and thus of a measurement 
of the absolute neutrino mass scale. However, in very re¬ 
cent years, interest in the double positron decay, /?■*"/?+, 
positron emitting electron capture, ECf3~^, and double 
electron capture ECEG, has been renewed. This is due 
to the fact that positron emitting processes have interest- 
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ing signatures that could be detected experimentally jH . 
In a previous article @ we initiated a systematic study 
of EGj3^, and ECEG processes and presented 

a calculation of phase space factors (PSF) for 2z//?"'"/?+, 
2vECj3^, 2i/EC EC and , 2i/i?C/?+. The process 

OvECEC cannot occur to the order of approximation 
used in since the emission of additional particles, 77 
or others, is needed to conserve energy and momentum. 
In this article, we focus on calculation of neutrinoless de¬ 
cay nuclear matrix elements (NME), which are common 
to all three modes, and half-life predictions for 0z^/?+/?"'" 
and 0vEC(3~^ modes. Results of our calculations are re¬ 
ported for nuclei listed in Table HI 


II. RESULTS 

A. Nuclear matrix elements 

The theory of Oz^/?/? decay was first formulated by 
Furry [l^ and further developed by Primakoff and Rosen 
[l5| , Molina and Pascual EH , Doi et al. Ea , Haxton and 
Stephenson EH > and, more recently, by Tomoda Ea and 
Simkovic et al. [23 ■ All these formulations often dif¬ 
fer by factors of 2, by the number of terms retained in 
the non-relativistic expansion of the current and by their 
contribution. In order to have a standard set of calcula¬ 
tions to be compared with the QRPA and the ISM, we 
adopt in this article the formulation of Simkovic et al. 
[^ . A detailed discussion of involved operators can also 
be found in Ref. jdj. 

We consider the decay of a nucleus into a nu¬ 
cleus example is shown in Fig. [T] If the 

decay proceeds through an s-wave, with two leptons in 
the final state, we cannot form an angular momentum 
greater than one. We therefore calculate, in this article, 
only Ozz/?/? matrix elements to final O"*" states, the ground 
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TABLE I. Double-/? decays considered in this article, the mass difference between neutral mother and daughter atoms, M{A, Z) — 
M{A, Z — 2), and their isotopic abundances. 


transition 

M(A, Z) - M{A, Z - 2)(keV)i^ 

P{%) 

ilNiao ^iFe32 

64Zn34 ^llNiae 

3gKr42 —>'34Se44 

I1Ru52 ^11Mo 54 
if Cdss Pdeo 

IfXejo ^lfTer2 

56°Ba74 —>-54° Xeyg 

gf Ceys —Bago 

1926.3 ±0.3 

1094.8 ±0.7 

2846.3 ±0.7 

2714.51 ±0.13^ 

2775.39 ±0.103 

2865.4 ±2.2 

2619 ±3 

2378.53 ±0.27^ 

68.077 ±0.009 
49.17 ±0.75 
0.355 ±0.003 
5.54 ±0.14 

1.25 ±0.06 
0.0952 ±0.0003 
0.106 ±0.001 
0.185 ±0.002 


Ref. 
Ref. 
<= Ref. 
Ref. 




FIG. 1. The decay foeCdsg ^loePdeo, an example of double- 
decay. 


state 0^, for which, in a previous article Q we have cal¬ 
culated the phase space factors, and to the first excited 
state 0^. 

In order to evaluate the matrix elements we make use 
of the microscopic interacting boson model (IBM-2) |2^. 
The method of evaluation is discussed in detail in Ref. [l| 
for double electron decay (/3“/3“). For double positron 
decay and positron emitting electron capture 

{ECP^) the same method ^plies except for the inter¬ 
change IT ^ 1 / in Eq. (5) of [l| and in the mapped boson 
operators of Eq. (18) of [Ij. The matrix elements of the 
mapped operators are evaluated with realistic wave func¬ 
tions, taken either from the literature, when available, or 
obtained from a fit to the observed energies and other 
properties {B{E2) values, quadrupole moments, B{M1) 
values, magnetic moments, etc.). The values of the pa¬ 
rameters used in the calculation are given in Appendix 
A. 

Here, we present our calculated NME for the decays of 


Table m The NMEs depend on many assumptions, in par¬ 
ticular on the treatment of the short-range correlations 
(SRC). In Table HI] we show the results of our calculation 
of the matrix elements to the ground state, 0)'', and to 
the first excited state, 0J, using the Miller-Spencer (MS) 
parametrization of SRC, and broken down into GT, F 
and T contributions and their sum as 

Mo. = 

We note that we have two classes of nuclei, those in 
which protons and neutrons occupy the same major shell 
{A = 64, 78,124,128,130,136) and those in which they 
occupy different major shells {A = 58,96,106). The mag¬ 
nitude of the Fermi matrix element, which is related to 
the overlap of the proton and neutron wave functions, 
is therefore different in these two classes of nuclei, be¬ 
ing large in the former and small in the latter case. This 
implies a considerable amount of isospin violation for nu¬ 
clei in the first class. This problem has been discussed 
in detail in Ref. Q and will form a subject of subsequent 
investigation. It is common to most calculations of NME 
and has been addressed recently within the framework of 
QRPAin Refs. Here we take it into account by 

assigning a large error to the calculation of the Fermi ma¬ 
trix elements. In the same Ref. Q it is also shown that 
the NME depend on the short range correlations (SRC), 
and that use of Argonne/CD-Bonn SRC increases the 
NME by a factor of 1.1-1.2. The same situation occurs 
for decay. In order to take into account the sen¬ 

sitivity of the calculation to parameter changes, model 
assumptions and operator assumptions , we list in Ta¬ 
ble [TTIl IBM-2 NMEs with an estimate of the error. The 
values of the Oj*" matrix elements vary between 2.3 — 6.1, 
the matrix element for the ®^Zn—>-®'^Ni transition being 
notably the largest. They are therefore of the same order 
of magnitude than the nuclear matrix elements for 

















3 


TABLE 11. IBM-2 nuclear matrix elements (dimensionless) for neutrinoless /ECP'^ /ECEC decay with Jastrow 

M-S SRC and qvIqa = 1/1.269. 


Nucleus 


0 ^- 




0 ^ 


M^ot 




M^ot 




58Ni 

2.072 

-0.152 

0.144 

2.310 

2.042 

-0.153 

0.101 

2.237 

®4Zn 

4.762 

-2.449 

-0.156 

6.127 

0.633 

-0.360 

-0.019 

0.837 

78Kr 

3.384 

-2.146 

-0.238 

4.478 

0.771 

-0.479 

-0.055 

1.014 

96Ru 

2.204 

-0.269 

0.112 

2.483 

0.036 

-0.012 

0.001 

0.045 

lo^Cd 

2.757 

-0.255 

0.191 

3.106 

1.395 

-0.110 

0.074 

1.537 

i24Xe 

3.967 

-2.224 

-0.192 

5.156 

0.647 

-0.359 

-0.032 

0.839 

i30Ba 

3.911 

-2.108 

-0.176 

5.043 

0.285 

-0.152 

-0.014 

0.366 

i36Ce 

3.815 

-2.007 

-0.161 

4.901 

0.318 

-0.167 

-0.014 

0.408 


decay, 2.0-5.4. 

In the same Table HIT] we also compare our results with 
the available QRPA calculations from Ref. [2^ with the 
addition of some more recent calculations from Refs. ^t\ . 
[H]. The QRPA [2^ NMEs are calculated taking into 
account GT and F contributions, and using the value 
gA = 1.25. As in the case of P~I3~ decay, QRPA tend 
to give larger values than IBM-2 and these two methods 
seem to be in a rather good correspondence with each 
other. 


TABLE III. IBM-2 matrix elements with M-S SRC and error 
estimate compared with available QRPA calculations. 


Decay 



oj 


IBM-2 

QRPA^ 

IBM-2 QRPA 

58Ni 

2.31(37) 

1.55 

2.24(36) 

64Zn 

6.13(116) 


0.84(16) 

78Kr 

4.48(85) 

4.19 

1.01(19) 

96Ru 

2.48(40) 

3.25 3.22-5.83^ 

0.05(1) 1.28-2.26^ 

lo^Cd 

3.11(50) 

4.12 5.94-9.08= 

1.54(25) 0.66-0.91= 

i24Xe 

5.16(98) 

4.78 

0.84(16) 

i30Ba 

5.04(96) 

4.98 

0.37(7) 

i36Ce 

4.90(93) 

3.09 

0.41(8) 


Reference 

Reference 

Reference 




B. Predicted half-lives for > 0)*' transitions 

The calculation of nuclear matrix elements in IBM-2 
can now be combined with the phase space factors cal¬ 
culated in jil to produce our final results for half-lives 
for light neutrino exchange in Table HVl and Fig. The 
half-lives are calculated using the formula 


where i = (3^(3^, EC. The values in Table HVl and 
Fig. [5] are for {rriu) = leV. They scale with {rriiPf' for 
other values. 


TABLE IV. Calculated half-lives in IBM-2 M-S SRC for neu¬ 
trinoless double-/?^ decay and positron emitting electron cap¬ 
ture for {ruv) = 1 eV and qa = 1.269. 


Nucleus 

/3+/3+ 

Ti/2 (102V) 

ECj3+ 

58Ni 


213 

64Zn 


52.9 

78Kr 

2.01 

0.79 

96Ru 

19.3 

1.70 

106cd 

10.8 

0.80 

i24Xe 

3.32 

0.19 

i30Ba 

15.4 

0.23 

i36Ce 

174 

0.27 



FIG. 2. (Color online) Expected half-lives for {niv) = 1 eV, 
QA = 1.269. The figure is in semilogarithmic scale. 
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Comparing the half-life predictions listed in Table HVl 
to the ones reported in Ref. Q for Oz//3“/3“ we can see 
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that values reported here are much larger. This is due 
to the fact that in cases studied here the available ki¬ 
netic energy is much smaller compared to decay. 

Furthermore, the Coulomb repulsion on positrons from 
the nucleus give s a smaller decay rate. As concluded also 
in Refs. (UllSSli the 0i^AC'/3“tdecay is expected to 

have the shortest half-live. In case of the neutrinoless 
double electron capture process, QvECEC, the available 
kinetic energy is larger and Coulomb repulsion does not 
play a role. However, this decay mode cannot occur to 
the order of approximation we are considering, since it 
must be accompanied by the emission of one or two par¬ 
ticles in order to conserve energy, momentum and angular 
momentum. 


III. CONCLUSIONS 

In this article we have presented evaluation of nuclear 
matrix elements in I3~^/OvECj3^/OvECEC within 

the framework of IBM-2 in the closure approximation. 
The closure approximation is expected to be good for 
these decays since the virtual neutrino momentum is of 
order 100 MeV/c and thus much larger than the scale of 
nuclear excitations. By using these matrix elements and 


the phase space factors of Ref. j^, we have calculated 
the expected 0vl3~^/0vECl3~^ half-lives in all nuclei of 
interest with gA = 1.269 and gy = 1, given in Table HVl 
and Fig. 
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IV. APPENDIX A 

A detailed description of the IBM-2 Hamiltonian is 
given in and [l^. For most nuclei, the Hamilto¬ 
nian parameters are taken from the literature fool-l^. 
The values of the Hamiltonian parameters, as well as the 
references from which they were taken, are given in Ta¬ 
ble El The quality of the description can be seen from 
these references and ranges from very good to excellent. 
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TABLE V. Hamiltonian parameters employed in the IBM-2 calculation of the final wave functions along with their references. 


Nucleus 


Ed, 


X!/ 


a 

a 

a 



(-l> 

C-TT 

cP 


SSj^ja 

1.454 














58pg a 

0.98 

0.98 

-0.26 

0.00 

-0.40 

0.80 

0.80 

0.80 







®4Zn [30| 

1.20 

1.20 

-0.22 

-0.25 

-0.75 

-0.18 

0.24 

-0.18 

-0.30 

-0.50 

0.30 

-0.30 

-0.50 

0.30 

a 

1.346 









-0.415 

0.082 





0.96 

0.96 

-0.18 

-0.495 

-1.127 

-0.10 


-0.10 







^®Se [32] 

0.99 

0.99 

-0.21 

0.71 

-0.90 



-0.10 







96Ru 

1.08 

1.08 

-0.21 

0.80 

0.40 

0.25 

0.25 

0.25 

0.30 

0.10 

-0.50 




®®Mo 

0.73 

1.10 

-0.09 

-1.20 

0.40 

-0.10 

0.10 

-0.10 

-0.50 

0.10 





106 Cd [3^ 

1.05 

1.05 

-0.325 

1.25 

0.00 

-0.18 

0.24 

-0.18 

0.20 

0.15 

0.00 




loepd [^ 

0.760 

0.844 

-0.160 

-0.22 

-0.30 

0.20 

0.05 

0.00 

-0.45 

-0.20 

0.01 




i24Xe 

0.70 

0.70 

-0.14 

0.00 

-0.80 

-0.18 

0.24 

-0.18 

0.05 

-0.16 





i24Te [34| 

0.82 

0.82 

-0.15 

0.00 

-1.20 

-0.18 

0.24 

-0.18 

0.10 






looBa 

0.70 

0.70 

-0.175 

0.32 

-0.90 

-0.18 

0.24 

-0.18 

0.26 






looxe 

0.76 

0.76 

-0.19 

0.50 

-0.80 

-0.18 

0.24 

-0.18 

0.30 

0.22 





looce [3^ 

0.90 

0.90 

-0.21 

0.79 

-1.00 

-0.18 

0.24 

-0.18 

0.26 

-0.11 





loofia 

1.03 

1.03 

-0.23 

1.00 

-0.90 

-0.18 

0.24 

-0.18 

0.30 

0.10 






^ Parameters fitted to reproduce the spectroscopic data of the 
low lying energy states. 







